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Abstract 

The top quark forward-backward asymmetry A^p^ measured at the Tevatron is above the Stan- 
dard Model prediction by more than 2(t deviation, which might be a harbinger for new physics. 

In this work we examine the contribution to A^p^ in two different new physics models: one is 
(N ■ 

' the minimal super symmetric model without R-parity (RPV-MSSM) which contributes to A^pn via 

!>■ ; 

■ sparticle-mediated t-channel process dd ^ tt; the other is the third- generation enhanced left-right 



X 



model (LR model) which contributes to A^p^ via Z'-mediated t-channel or s-channel processes. 
We find that in the parameter space allowed by the tt production rate and the tt invariant mass 



' distribution at the Tevatron, the LR model can enhance Ap^ to within the 2a region of the Teva- 



tron data for the major part of the parameter space, and in optimal case Ap^ can reach 12% 



^ I which is slightly below the la lower bound. For the RPV-MSSM, only in a narrow part of the 

parameter space can the A" couplings enhance A^p^ to within the 2a region while the A' couplings 
just produce negative contributions to worsen the fit. 

PACS numbers: 14.65.Ha,14.80.Ly,11.30.Hv 
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I. INTRODUCTION 



As the heaviest fermion with a mass at weak scale, topquark is speculated to be a sensitive 
probe of new physics beyond the Standard Model (SM) The precision measurement of its 
properties, now being performed at the Tevatron and will be continued at the CERN Large 
Hadron Collider (LHC), will either unravel or further constrain the new physics related to 
the top quark. 

So far the production rates of top pair and single top measured at the Tevatron are in 
good agreement with the SM predictions, but a more than 2a deviation is reported in the 
forward-backward asymmetry A^p^ in top pair production, which is defined by 

A^t(cos^ > 0) -iVi(cose < 0) 



APP 



(1) 



Nt{cose > 0) + Nt{cose < 0) ' 

with 9 being the angle between the reconstructed top quark momentum and the proton beam 
direction in ti rest frame. In the SM, this asymmetry gets dominant contribution from the 
next-to-leading-order QCD correction and was found to be several percent: Ay-q{^M) = 
5.0 ± 1.5% Compared with its experimental value A^-Q{exp) = 19.3 ± 6.9% measured by 
the CDF and DO collaborations , this SM prediction is below the measured value by more 
than 2a deviation. 

Such a discrepancy might be a new physics footprint in top quark sector and has been 
studied in several new physics models, where the Kaluza-Klein excitations in extra dimen- 

n p ^ ' 

sions 1^, the presence of new gauge bosons {Z', W', axigluon) [5|, |6[ or new scalars [7| are 
utilized to try to explain the discrepancy. Noting that the mechanisms proposed in these 
literatures can also be realized in some popular new physics models, we in this work study 
the asymmetry A'p^ in the supersymmetric models and the left-right models jsj. 

For the minimal supersymmetric standard model (MSSM) [9], the SUSY influence on tt 
production comes from loop effects, in which the SUSY-QCD effects are dominant over the 



SUSY-EW effects 



10|. But among the SUSY-QCD one-loop diagrams only the box diagrams 



contribute to the asymmetry A^p^ and, consequently, the contribution is negligibly small (see 
the discussion in ^). Therefore, in our analysis we consider the R-parity violating MSSM 
(RPV-MSSM) ll| which allows for tree-level contribution from t-channel process dd — )■ ti by 
exchanging a color-singlet slepton or a color-triplet squark. For the general left-right models, 
since the predicted new gauge bosons are usually at TeV scale and unlikely to affect Ap^ 
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significantly, we here consider a special left-right model called the third-generation enhanced 
left-right model jl2|. This model predicts a new gauge boson Z' which contributes to 
via t-channel or s-channel processes. 

This paper is organized as follows. In Sec. II we describe the calculation of the asymmetry 
A^pj^ in the RPV-MSSM and present some numerical results and discussions. In Sec. Ill 
we perform similar analysis in the third-generation enhanced left-right model. Finally, some 
discussions and the conclusion are presented in Sec. IV. 



II. A^p^ IN i^-PARITY VIOLATING MSSM 

In the popular MSSM, the invariance of -R-parity, defined by i? = (_ i)25+3B+l ^ field 
with spin S, baryon-number B and lepton-number L, is often imposed on the Lagrangian 
in order to maintain the separate conservation of B and L. Although i?-parity plays a 
beautiful role in the phenomenology of the MSSM (e.g., forbid proton decay and ensure a 
perfect candidate for cosmic dark matter), it is, however, not dictated by any fundamental 
principle such as gauge invariance and there is no compelling theoretical motivation for it. 
The most general superpotential of the MSSM consistent with the SM gauge symmetry and 
supersymmetry contains i?-violating interactions which are given by pJJ] 

^9 = 2^ijkLiLjE^ + X'ijf^LiQjDl + -A-'j-^e^^'^L^^^Dj^D^^ + fiiLiH2, (2) 

where k are generation indices, c denotes charge conjugation, a, f3 and 7 are the color 
indices with e"^"^ being the total antisymmetric tensor, H2 is the Higgs-doublet chiral super- 
field, and Li{Qi) and Ei{Ui, Di) are the left-handed lepton (quark) doublet and right-handed 
lepton (quark) singlet chiral superfields. The dimensionless coefficients Xijk (antisymmetric 
in i and j) and in the superpotential are L- violating couplings, while A-'^-^ (antisymmetric 
in j and k) are 5-violating couplings. 

The expression of implies that both A^^-^^and can induce new top quark interac- 
tions. In terms of the four- component Dirac notation, the interactions involved in top pair 
production are 

= Kj'iA< - iKS'AdL + d^Ad'n + h.c, (3) 
and these interactions can contribute to the forward-backward asymmetry by the diagrams 
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FIG. 1: Feynman diagrams contributing to A^p^ in the RPV-MSSM, where and dj^ denotes a 
left-handed slepton in i-th generation and right-handed squark in A;-th generation, respectively. 

TABLE I: The upper bounds on the couplings A^g^ (i = 1,2,3) and Ag'^^ {k = 2,3) 



couplings 




bounds 


sources 




0.03 


m-.^/(100 GeV) 


QwiCs) 


^^231 


0.18 


mjfe/(100 GeV) 




\' 

^^331 


0.26 


mjfc /(lOO GeV) 


K —7- TTVU 


a" 

^31fc 


0.97 


mjfc /(lOO GeV) 




a" 

^31fc 




1.25 


perturbativity 



shown in FiglH The corresponding amphtudes are then given by 

mRpv\ ,-a a I\' \2 u{t)PRu{d)v{d)PLv{t) 

n/rRPV\ _ V . i\" \2 u{t)-1^PBv{t)v{d)-1^'PRu{d) 

^dd^m" - -^^/JpA^aaAlAgifcl 2[(j9i -^4)2 - ] ' 

with a, /3, p, cr and A being color indices of the quarks and squarks. These amphtudes affects 
A^p^ by interfering with the QCD amphtude dd ^ g* ti and also by its own square. In 
our results presented below, we have included the SM contribution to A^p^ and considered 
only one coupling non-zero each time. 

The SUSY parameters involved in the calculation are the couplings A^g^^ and Ag^^^ as well 
as sparticle masses. So far both theorists and experimentalists have intensively studied the 



13| and obtained some bounds 



Q. 



phenomenology of these couplings in various processes 
In Table [I] we list the relevant bounds and, as can be seen, in case of heavy squarks, these 
bounds are quite weak. Note that for Ag^;^ we do not use the stringent bound from n — n 



oscillation [ij] because they depend on additional SUSY parameters which are not involved 
in our processes. 

In Figj2] we show the dependence of A*p^ on the relevant sparticle masses. Here we 
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FIG. 2: A^p^ versus the sparticle mass in the R-parity violating MSSM. The left (right) panel is 
induced by X[^^ (-^s'lfc) which are fixed at their maximally allowed values. The SM contribution to 
A*p^ is also included. 

assume all squarks and sleptons are degenerate and sum over the contributions from different 
generations. The couplings A^g^ and Ag^^ are fixed at their maximally allowed values which, 
as shown in Table I, vary with squark mass. The SM parameters are taken as ISj 



mt 



172.5 GeV, mz = 91.19 GeV, sinOw = 0.2228, a,(mt) = 0.1095, a = 1/128. (6) 



Figj2] indicates that both X[^^ and Ag^^ can give a negative contribution to Ap^ and in 
the worst cases, they decrease A^^^ by 2% and 5% respectively, which are comparable in size 
with A^p^^^ . In a narrow part of the squark mass (250 ~ 400 GeV), the couphng Ag^^ can 
also give a positive contribution to enhance A'^p^ to within the 2a region of the Tevatron 
data. 

We note that the scalar-mediated contributions have been analyzed in a model indepen- 
dent way in 3] and the results were presented for color-singlet, -triplet, -sextet and -octet 
scalar respectively. We checked that our analytic results are in agreement with [3] for the 
color-singlet and -triplet cases except that our study is restricted in a specified model, the 
RPV-MSSM. 

Before en ding this section, we give a comment on A^^^ in the top-color assisted technicolor 
model (TC2)[16]. This model predicts some composite bosons, 7r° and nf, with mass at weak 
scale and having large Yukawa couplings to top quark. As a result, Ap^ gets additional 
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contribution from t-channel processes uu{cc) — > tt by exchanging a color-singlet scalar vr^, 
which is similar to Fig{T](a) in RPV-MSSM. Noting the up quark content in proton is larger 
than the down quark content, one may expect a larger effect on A^p^ in TC2 model than 
in RPV-MSSM. This is not true because for uu — > tt in TC2 model, its contribution is 
proportional to ur — t^ mixing which is determined by the triangular texture of the up-type 
qua* m^. matrix and tu™= out to be very small R . As to ^ ti, although mWng 
may be sizable [iTi], it actually gives no contribution to A^p^ because the distributions of c 
and c in proton are approximately same and so the initial state of this process is symmetric 
under the exchange of c and c. In fact, we numerically calculated all TC2 contributions 
to tt production at the Tevatron, which include the s-channel processes gg,bb 7r°* — > tt, 
t-channel processes bb — > tt by exchanging and uu, cc — t- tt by exchanging vr^, and we 
found that TC2 model can change au by at most 200 fb and A'^p^ by order of 10~^. 



III. A^pg IN THE THIRD-GENERATION ENHANCED LEFT-RIGHT MODEL 



In the third-generation enhanced left-right model 12|], the gauge group is SU{3)c x 
SU{2)l X SU{2)ii X U{1)b-l with gauge couplings g^, gi, gn and g, respectively. This 
model differs from other left-right models by having right-handed gauge bosons couple pre- 
dominantly to the third generation fermions. Noting that the right-handed gauge boson Zr 
will mix with the standard Zq to form mass eigenstates Z and Z', one can write down the 
neutral gauge interactions of quarks as 



I cos Piy ^ 



where tan6'/j = gjgn, gy = gcosOji = gRsinOn, is the mixing angle between Zr and Zq, 
and V^^ are unitary matrices which rotate the right-handed quarks u^n and dju from inter- 
action basis to mass eigenstates. This model also predicts new charged gauge interactions 
of quarks, but since their effect on A^^^ is small, we do not consider them here. Note that 
in Eq.(I71), q and ql are summed over all quarks, and the repeated generation indices i and 
j are also summed. 

Eq-dZ]) indicates that the Z'uiUj interaction is large when gji ^ gy or cot 9 ji ^ 1. 



U{pi) tips) U{pi) t{p3) u{pi) t{pz) U{pi) t{p3 
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FIG. 3: Feynman diagrams contributing to A^p^ in the third-generation enhanced left-right model. 

This feature can be utilized to enhance A^pg by the t-channel process uu — )■ ti if {y^)ut is 
moderately large. In [12I, and are assumed to be nearly diagonal to satisfy the severe 
flavor- changing neutral-current (FCNC) constraints. However, as pointed out in [6|, a sizable 
UR—tfi mixing does not conflict with these constraints given that the other flavor mixings are 
suppressed. Here we further point out that this pattern of flavor mixings does not necessarily 
require the up-top element in up-type quark mass matrix M„ to be much larger than other 
off-diagonal elements. For example, assuming (V^)^* = 0.2, (V^)^ = and (V^)mc = 0, we 
numerically solve the equation M^MuV^ = Mj^^^ with Mj-^^ = diag{ml^,ml,m'f}, and 
we find it possible that {Mu)ct is several times larger than {Mu)ut- 

In the third-generation enhanced left-right model, beside the dominant QCD contribution 
qq ^ g* ^ tt, diagrams shown in Figj3] also contribute to ti production at the Tevatron, 
and the amplitudes of these new contributions are 

, f ^ V ^(^)7A.[gk-PL + ffkPi;]t^(t)t;(n)7^[gkPL + g|fiPi;]n(tx) 



M, = t5^p5,. —— -'--^ , Y 2 -r " (9) 

A/f -A X ( ^ ( I ^- a m2|t/« t/« ^2 Kt)ltiPRu{u)v{u)-^>'PRv{t) 

Mc = idapdpaii: [(,zSnj[coteR + tanen)] \VjituVRtt\ 7 2 '(^0) 

\ZCu]S^ / {^pl p^j 



where Tz' is the Z' width obtained by adding all decay channels of Z', = sin ft 
cosOy^, and the coupling coefficients 5'^'^, (7^^, (7^"^ and (7^"^ are defined as 

9zL = ^ - - ^s^ tan Or^z, (12) 
4 4 

9zR = -^sl- -s^tanOR^z, (13) 
4 1 

9zR = ~3 ~ s'^'" ^^^^ + (-^^) 
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^z^L = i'^ - ^^D^z + ^s^ tan 9r, 

4 4 
9z'R = -^sl^z + ^s^tanOR, 



4 1 

■ 3 s^^z + tan 9r^z - Su> cot 6'ij. 



(15) 
(16) 
(17) 



Note that among the four amphtudes, only Mc and Md interfere with the QCD amphtude. 

In this model the new parameters cot Or, M^i and (V^)„.f are involved in our cal- 
culation. Constraints on these parameters were discussed in 12(|, and it was found that 
< < 0.02 and cot Or < 20. As for Mz', we should note that the constraints from 
CDF search for Z' ISj and from the global fitting of the electroweak precision data [19i] are 
invalid here since these constraints arise mostly from the processes involving the first- or 
second-generation of fermions. So far the pertinent bound comes from e~^e^ — ?■ bb at LEP-II, 
which requires Mz' ^ 460 GeV for cot Or > 10 12|. In following analysis, without stating 
explicitly, we include the QCD contribution to A^p^ and scan the new free parameters in the 
following ranges 



500 GeV < Mz' < 2000 GeV, < < 0.02, 10 < cot Or < 20, 0.1 < {V^)ut < 0.2 



- s . 




Apg''= 19.3 ± 6.9 % 

a: 
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FIG. 4: Scatter plots of A^p^ versus cot Or and m^/ for the no-mixing case in the LR model. 



In our discussion we consider two cases, i.e., with and without ur — Ir mixing. For the 
no-mixing case, new contributions to A^p^ only come from diagrams (a) and (b) of Figj3] 
and the dependence of Ap^ on cot Or and m^' are shown in FigJH This figure shows that a 
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light Z' with relatively small cot Qb. can enhance J^p^ ^^e 2a region of its experimental 
value. This is because small cot can enhance Z'uu interactions, and a light Z' can make 
the resonance effect of the diagram Figl2]^b) on tt production more significant. FigJH also 
shows that in the no-mixing case, the Z' contribution to J^p^ can only reach 1%. This is 
due to the smallness of Z'uu couplings (see Eg. (1141) and Eq. (fT5l) ) and the absence of the 
interference of diagram (a) and (b) in Figj3]with the dominant QCD amplitude. 
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FIG. 5: Scatter plots of Ap^ versus cot 9r in the mixing case of the LR model. The left (right) 
panel is without (with) considering the constraints from and M^f. 




In Figsl5]and[6l we show Ap^ versus cot 9r and Mz' respectively in the mixing case. The 
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left panels of these figures show that the effect of Figj3] (mainly diagram (d)) is potentially 
very large, pushing the value of A^p^ up to 50%. In this case, the effects on the total ti 
production rate and the tt invariant mass distribution are also large and so it is necessary 
to consider constraints from these observables measured at the Tevatron. Now based on the 
CDF 4.6 fb~^ luminosity data, the measured total cross section is 0"°^^ = 7.50 ± O.Slgtat =t 



O.SAsyst ± 0.15f/i pb for nit = 172.5 GeV 20|]. Combining errors in quadrature, one can get 



'^a^ = 7.50 ± 0.48 pb, which is in good agreement with the SM prediction au = 7-5~^o^ 



pb j2l|. The invariant mass distribution was also measured by CDF, and the results are 
presented in nine bins of Mu When we calculate these observables, we multiply an 

overall K-factor 1.329 for the dominant tree-level QCD contribution ^l|. The right panels 
of FigsIS] and IH] are then obtained by requiring the total cross section and the differential 
cross section in each bin to be within the 2a regions of their experimental values. These 
results shows that even with the constraints, the third-generation enhanced left-right model 
can still enhance A^p^ to 12% (well above the 2cr lower bound), but it can not enhance A^^^ 
to within the la region of its experimental value. Note that unlike the no-mixing case, large 
contribution to A^p^ comes from the region where cot 6^ is large. The reason is that in the 
mixing case, dominant contribution arises from diagram (d) of FigJSl and this contribution 
is proportional to (cot 9r + tan 6'/?)^. 

IV. DISCUSSION AND CONCLUSION 



From our study and the previous works {2!, Q-ll, |23|], we can learn that to make sizable 
contribution to A^p^, the following two conditions must be satisfied. One is the initial state 
must be uu and/or dd and in case that only dd initiated contribution is responsible to explain 
A^p^, the involved interaction must be strong enough to compensate the suppression of the 
parton distribution of down quark in proton. The other is the amplitude of the ti production 
must contain terms proportional to cos 6 with 6 being the angle between the reconstructed 
top quark momentum and the proton beam direction in ti rest frame, or in other words, 
contain terms like {pu ■ Pt){Pu ■ Pt) — {Pu ■ Pt){Pu ■ Pt)- This requirement implies that new 
physics affect ti production through the following ways: 

(1) Through s-channel process qq — )■ ti by exchanging a gauge boson 0, S]. If this process 
interferes with the the s-channel QCD process gg — t- (7* — )■ ti, the interaction of 
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the gauge boson with q and t must have axi- vectorial component. Examples in this 
direction are the presence of Kaluza-Klein excitation of gluon in extra dimension jj] or 
the axigluon ^. If this process does not interfere with the s-channel process, to affect 
A^p^ by itself, the interaction must have both vectorial and axi- vectorial component, 
like what we studied in diagram (b) of Fig|31 From yet known studies we can infer 
that it seems difficult for the latter case to enhance A^p^ significantly without spoiling 
the constraints from ati and Mti. 

(2) Through t-channel process qq — )■ tt by exchanging a vector boson or a scalar 

a a, 

which interferes with the QCD process qq ^ g* ^ tt. In this way, scalar is less 
efficiency than vector boson in explaining A^p^ given that they have the same coupling 
strength and mass. This is because for the scalar case, there is a competition between 
spin correlation and the Rutherfold singularity Moreover, as shown in fzl, the 
scalar-mediated contributions can be categorized by the transformation property of 
the scalar under SU{3)c- If the scalar is color-triplet or -sextet, there exists large 
parameter region to explain A^p^ within la and at the same time to remain au within 
the experimental errors, while for color-singlet or -octet scalar, it is very difficult to 
jroduce a large positive contribution to A^p^ without spoiling the constraint from au 
7| . In our work, we checked this conclusion for color-singlet and -triplet cases. 

From our study, we can also learn that, although the top quark pair productions at the 
Tevatron and the LHC may be sensitive to new physics, the effects of new physics (like the 
popular MSSM or TC2 models) are usually not so large to be well above the experimental and 
theoretical uncertainties [l^. A complementary or even more sensitive probe for new physics 
effects in top quark sector is top quark FCNC processes, which are extremely suppressed and 
unaccessible in the SM but can be greatly enhance d by several orders to reach the observable 



24] or the TC2 model [25 1. 



level in some new physics models like the MSSM 

In summary, we in this work calculated the new physics contribution to the top quark 
forward-backward asymmetry A^p^ at the Tevatron in two different models: the minimal 
supersymmetric model without R-parity (RPV-MSSM) and the third-generation enhanced 
left-right model (LR model). We found that in the parameter space allowed by the ti produc- 
tion rate and the tt invariant mass distribution at the Tevatron, the LR model can enhance 
A^p^ to within the 2a region of the Tevatron data for the major part of the parameter space. 
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and in optimal case A^p^ can reach 12% which is shghtly below the la lower bound. For the 
RPV-MSSM, only in a narrow part of the parameter space can the A" couplings enhance 
Ap^ to within the 2a region while the A' couplings just produce negative contributions to 
worsen the fit. Noting that the R-parity conserving interactions in the MSSM are unlikely to 
give large enough contribution to Ap^, we conclude that the MSSM with (without) R-parity 
will be disfavored (favored in case of A" couplings) if the discrepancy of A^p^ persists with 
more forthcoming data accumulated by the Tevatron. We also checked the top-color assisted 
technicolor model and found that it gives negligibly small contributions to A^^^ and thus 
unlikely to explain the Tevatron data. 
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